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DEFINITION  OF  SYMBOLS  USED 
A = Absorption  Coefficient  in  dB/kyd. 

B = Leakage  Coefficient  in  dB/kyd. 

f = Frequency  in  kHz . 

H =•  Duct  depth  in  feet, 

n = Sea  State. 

R = Ran.-e  in  kiloyards. 

t = Temperature  (decrees  C). 

T = Temperature  (decrees  F). 
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1. 


A NEW  FORMULA  FOR  CALCULATING  ACOUSTIC  PROPAGATION  LOSS 

IN  A SURFACE  DUCT  ' 

preTcis 


1.  This  note  describes  the  development  of  a formula  for  the  calculation  of  the 
acoustic  propagation  loss  in  a surface  isothermal  duct.  The  formula  derived  is 
an  empirical  fit  to  a number  of  experimental  measurements  which  embrace  a limited 
range  of  the  factors  affecting  acoustic  propagation  loss,  and  extrapolation  out- 
side these  limits  should  obviously  be  regarded  with  caution.  The  variation  of 
propagation  loss  with  depth  within  the  duct  haB  not  been  investigated  due  to  lack 
of  data. 

2.  Comparison  of  the  measured  results  with  the  new  formula  indicate  a mean 
error  of  typically  less  than  t 2 dB,  with  a standard  deviation  of  7 dB  on  indivi- 
dual readings.  A total  of  438  measured  propagation  loss  values  has  been  used 
from  four  different  sources. 


UNLIMITED 


3. 


INTRODUCTION 


It  had  teen  found  ti.at  the  latest  version  of  the  formula  given  by  Saxton 
(Ref.  l)  gave  duct  propa^t  ion  losses  which  were  greater  than  the  measured 
values  on  a number  of  occasions.  For  this  reason,  a systematic  study  of  the 
differences  between  measured  values  and  those  predicted  by  the  Saxton  formula  was 
undertaken. 


PROPAGATION  LAW 

Using  the  symbols  already  defined,  the  simolified  Saxton  duct  propagation 
loss  mav  be  written:- 

Loss  (dB)  ~ 20  lor  R + 60  + (A  + B)  R 

O.c; 

for  short  ranges  (R  5 (0.122  H)  ) 
or  Loss  (dB)  = 10  lor  R + 5 log  H + 50.9  + (A  + B)  R 

0 5 

for  long  ranges  (R  > (0.122  H) 

5.  This  note  is  concerned  with  the  values  of  the  coefficients  A and  B 
(Absorption  and  leakage  coefficients). 

ABSORPTION  COEFFICIENT 


6.  The  coefficient  used  by  Sacton  was:- 
A = 2#  r (dB/kyd) 


7.  In  the  present  formula,  a more  recent  expression  derived  by  Hall  (Ref.  2) 
has  been  used:- 

.2 


1.776  f 


1.5 


1 


IITlM 


32.768  + f3  1 + ^4 

f 


0.65053  ft  f 0.026847  f" 

r— + — : — 


rt2 + r 


where  = 21.9  x 10 


( 6t  + 118^ 
\t  + 273  > 


8.  Fir.  1 shows  a comparison  of  these  two  absorption  coefficients  over  the 
frequency  range  1 to  10  kHz  for  temperatures  of  40°F  and  60°F.  The  new  formula 
is  cumbersome  to  use,  but  does  fit  experimental  data  much  more  closely  than 
previous  formulae. 


LEAKAGE  COEFFICIENT 


9.  Saxton’s  expression  for  this  coefficient  has  been  given  in  various  forms, 
a commonly  used  one  being 

B = K l0S  , (l.4)n  (dB/kyd) 

173  (H)°*^ 

where  K is  an  empirical  constant  whose  value  has  been  changed  from  time  to  time, 
a typical  value  being  1600. 
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4. 


UNLIMITED 


10.  Usinr  this  expression  for  the  leakage  coefficient  and  the  new  absorption 
coefficient,  theoretical  duct  propagation  losses  were  compared  with  four  sets  of 
experimental  data,  comprisinr  a total  of  438  measurements.  The  main  characteris- 
tics of  the  data  (ranre  of  parameters  encountered)  are  /riven  in  Appendix  I. 

11.  For  convenience  of  presentation,  the  differences  between  theoretical  and 
measured  loss  values  have  been  grouped  into  quasi-loparithmic  ranpe  classes,  and 
the  mean  error  has  been  dotted  at  the  mean  ranee  of  the  class  (Fir.  2).  'Error 
is  here  defined  as  the  measured  loss  subtracted  from  the  calculated  loss. 

1 .j  jt  be  note  1 that  there  is  an  obvious  ranine-dependent  component  in  the 

mean  errors  for  all  four  sets  of  data,  and  that  additionally  the  errors  for  the 
7.c,  data  are  nirher  than  those  for  3.2R  kHz  and  3.5  kHz. 

13.  Tr.e  ran/re-dp  pendent  error  component  indicated  that  the  empirical  scalin* 
factor  in  the  coefficient  B (=  1 -00 ) required  modification  to  obtain  a rood  fit 
to  the  data;  the  frequency -dependent  error  component  sufrested  that  the  "lor  f 
factor  in  the  coefficient  should  be  chanred.  A linear  frequency  dependence 
seemed  appropriate,  since  this  choice  also  removes  the  anomalous  effect  if  the 
formula  is  used  below  1 kHz. 

14.  The  denominator  of  the  coefficient  B was  also  chanred  to  ta/.e  into  account 
the  r-lirht  variation  of  ray  curvature  due  to  temperature,  as  in  the  ori/rinal 
Saxton  expression,  but  this  had  little  effect  compared  with  the  other  changes 
detailed  above. 

1r.  As  a result  of  these  three  changes  the  final  expression  chosen  for  the 
leakare  coefficient  is:- 

B = &LL — (l.4)n  (dB/kyd) 

((1452  + 3.5t)  H)U* 

16.  Fir.  3 shows  the  mean  errors  obtained  usinp  the  new  formula, and  it  is 
apparent  that  the  formula  rives  a much  better  fit  to  the  experimental  data.  It 
should  be  emphasised  here  that  the  frequencies  used  were  3.25,  3.5  and  7.5  kHz, 
and  that  the  ranre  of  environmental  conditions  is  limited  as  detailed  in 
Appendix  I,  extrapolation  outside  the  ranre  of  conditions  considered  may  invali- 
date the  formula,  especially,  for  instance,  vrith  a shallow  duct  at  low  frequencies. 
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APPENDIX  I:  SUMMARY  OF  EXPERIMENTAL  DATA 


Data  1 

Data  2 

Data  3 

Data  4 

Frequency  (kHz) 

3.25 

3.5 

3.5 

7.5 

Number  of  Readings 

174 

15 

159 

90 

Range  of  Duct  Depth  (feet) 

80  to  220 

130 

140  to  200 

80  to  210 

Range  of  Sea  State 

2-3  to  4-5 

2 

3 

2-3  to  4-5 

Ran re  of  Temperature  °F 

64  to  76 

64 

59  to  62 

65  to  76 

Ranre  of  Ranges  (kyds) 

1 to  34 

O 

4 to  30 

1 to  28 
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ABSORPTION  COEFFICIENT  <JB/  kyd 


MEAN  ERROR  <1B 


UNCLASSIFIED  UNLIMITED 


FIGS  2 & 3 


FIG  2 MEAN  ERRORS  USING  SAXTON  EQUATION 


-X  3-25  kH  z DATA  I 
35kHz  DATA  2 
+ 3-5  kH  z DATA  3 
-O  7-5  kHz  DATA  4 
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FIG. 3.  MEAN  ERRORS  USING  NEW  FORMULA 
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